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Abstract 

We study heat-bath Glauber dynamics for the ferromagnetic Ising 
model on a finite cycle (a graph where every vertex has degree two). We 
prove that the relaxation time T2 is an increasing function of any of the 
couplings J xy . We also prove some further inequalities, and obtain exact 
asymptotics for T2 at low temperatures. 

1 Introduction 

Consider the ferromagnetic Ising model on a finite graph G = (V,E). Let n 
denote the number of vertices. We allow variable interaction strength, so each 
edge e = [x, y] G E has a coupling constant J xy > 0. As usual, every vertex 
x G V has a spin a x G {+1,-1} and the probability of a configuration is given by 
the Gibbs equilibrium measure 7r(er) = Z^ 1 exp(J^ a . y \ eE Jxy^x^y)- We assume 
no external field and no boundary conditions. 

We are interested in the Glauber dynamics for this model. The Glauber 
dynamics is a Markov chain on the space of configurations of the Ising model, 
whose stationary distribution is the Gibbs measure. It comes in several flavors 
(heat-bath, Metropolis, etc.), but they follow the same scheme. At each step 
of the algorithm, one vertex is picked according to some rule (e.g. uniformly 
random) and with a certain probability one flips its spin (from +1 to —1 and 
viceversa) . The flip probability depends on the value of the spin and its neigh- 
bors; if it is chosen appropiately, then the chain distribution converges to the 
Gibbs measure. 

For any configuration a we let a x be the configuration obtained by flipping 
spin x, while leaving all other spins unchanged. We shall focus on the heat-bath 
Glauber dynamics, in which a vertex x is picked uniformly at random, then 
flipped with probability ^(a x ) / '(tt((j) + tt((t x )). Formally, this is the discrete- 
time Markov chain whose transition matrix A has size 2™ x 2™ and is defined 
by: 

A(a,a x ) = (l/n)ir(a x )/(-K(a)+Tr(a x )) 
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A(a,r) 



for all other t 
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One fundamental quantity is the relaxation time t 2 , denned by 



t 2 = (n(l - /i 2 )) 
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where ^t 2 is the second largest eigenvalue of A. We scale by a factor of n so 
that the r 2 defined above is the same as for a continuous-time chain where each 
vertex is updated at rate one. Essentially, it says we do n moves of the discrete 
chain per unit of time. 

In this paper we study the case when the graph G consists of only one cycle 
(every vertex has degree two) . Our main result is 

Theorem 1 For the ferromagnetic Ising model on the cycle, the relaxation time 
r 2 is an increasing function of any of the couplings J xy . 

This result was motivated by trying to understand how the dynamics de- 
pends on the parameters of the model. The Glauber dynamics has been used in 
computer simulations of the Ising model, and in describing the non-equilibrium 
behavior of the model. There has been a fair amount of research on its relax- 
ation time and mixing time; see |S] for a comprehensive survey, and |3] for an 
analysis of the Glauber dynamics on trees. Most existing results give asymp- 
totic estimates for the relaxation time when the graph size is large; in particular, 
there are fairly sharp estimates for trees and for the lattice Z d . However, much 
less is known about how the relaxation time depends on the parameters of the 
model, such as the temperature or the coupling constants. Such results would 
allow comparing mixing times for different graphs at different temperatures. 

Computer simulations and heuristics suggest that, at least in the ferromag- 
netic case, such dependence should be monotone. Indeed, consider the dynamics 
started at an initial state where all spins are +1. To achieve mixing, lots of spins 
would have to flip to —1. But in the ferromagnetic model, neighboring spins 
tend to have the same value. The higher the interaction strength, the harder it 
will be for spins to flip. This led to the following 

Conjecture 2 (Yuval Peres, personal communication.) For the ferromagnetic 
Ising model on any graph, the relaxation time r 2 is an increasing function of 
any of the couplings J xy . 

What this says is that for two sets of couplings J, J' with J xy < J xy , we 
must have t 2 (J) < t 2 (J'). At first glance it may seem that one way to prove 
this is to find a coupling (not to be confused with the couplings J xy ) between 
the two Glauber dynamics, in which the dynamics for J is dominated by the 
one for J'. This would imply the desired inequality for relaxation times. 

However, such a coupling does not exist. To see why, choose some starting 
state (for example, all spins are +1) and let P n ,P' n be the measures obtained 
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after n steps of the dynamics. Let n, n' be the corresponding Gibbs measures. 
We would like to define both dynamics on the same probability space so that 
stochastic dominance P n -< P' n holds. But then it must also hold in the limit, 
so we would need to have 7r -< 7r', which cannot hold, because both measures 
have mean zero. 

In this paper, we prove the conjecture in the special case when G is a cycle. 
The proof proceeds as follows. We show that on the cycle, the transition matrix 
of the Glauber dynamics preserves linear functions of the spins, and furthermore 
that there is a linear eigenfunction corresponding to the second largest eigen- 
value. Hence by looking at the action of the dynamics on the space of linear 
functions of spins, the original 2™ x 2™ transition matrix can be reduced to a 
smaller, n x n matrix, whose entries are easily computable, and mostly zeroes. 

We analyze the smaller matrix to obtain a variational description of the 
relaxation time, as well as an inequality involving the components of its lead 
eigenfunction. These two combine to prove the desired monotonicity. In addi- 
tion, we obtain some further inequalities involving the relaxation time, as well 
as low-temperature asymptotics. 

2 Two Useful Lemmas 

The results in this section are valid for the ferromagnetic Ising model on any 
graph. They will be needed later in the proof of our main theorem. 

The transition matrix A of the heat-bath Glauber dynamics was defined 
above. It acts on the space of functions of u in the natural way. It has 2™ 
eigenvalues 1 = ii\ > fj,2 > • • • > M2" with corresponding eigenfunctions f-y = 
1, / 2 , . . . , /2". Because A is 7r-reversible we have ir(fifj) = for i ^ j (if an 
eigenvalue is multiple, we can choose its eigenfunctions to be 7r-orthogonal) , and 
in particular tt/, = for i > 2 (we use the standard notation nf = E^f). 

An important class of functions are the increasing functions, defined in the 
usual way: / is increasing if f(a) < /(r) whenever a < t (that is Oj < r^Vi). / 
is strictly increasing if strict inequality holds. We have 

Lemma 3 (Y. Peres) The second eigenvalue /12 has an increasing eigenfunc- 
tion. 

Proof. The method is similar to the proof of Perron- Frobenius. We need 
a standard result about the Glauber dynamics: if / is increasing then Af is 
increasing (for a proof see, for example, Theorem 2.14 in 0]). 

Let / be any increasing function with wf = 0; it can be written as a linear 
combination of eigenvectors / = Y)i—o Qifi- Let // = max(|/i 2 |, |/Z2« |). Iterating, 
we have 

2™ 

A m f/[i m = J2^/ri m fi ( 3 ) 

i=2 
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so if we let 

2" 

9 = ^QifdifJ-i = M) 

i=2 
2™ 

i=2 

we obtain convergence A 2m f / ' \j? m — > <? + h, A 2m+1 f / fi 2m+1 — > 5 — ft, where 
51 is some eigenvector corresponding to /Z2 (a multiple of /a if /12 is a simple 
eigenvalue). Since A n f is increasing for all n, the limits must also be increasing, 
so their sum must be increasing, so g must be increasing. 

It remains to show that g is nonzero. If it were then we would have both h 
and —h increasing, so h would also have to be zero. To avoid this, it suffices to 
choose / not orthogonal to /a and /an. For example, / = /a + h n + K X}* 7 * 
has 7r/ = 0, and for large if > 0, 7r(//a) 7= 0, 7r(//a») / 0, and / is increasing. 
Hence 92 7^ 0, (72" 7^ so at least one of <? and /i is nonzero. 

Lemma 4 fy. Peres) If A has a strictly increasing eigenfunction f , then f 
corresponds to fi2- 

Proof. Let g be any other increasing eigenfunction. Wc will show there are 
only two possibilities: g is constant (so its eigenvalue is 1), or g corresponds to 
the same eigenvalue as /. Otherwise, we would have ir(fg) = 0, nf = 0, ixg = 0. 
Since / is strictly increasing, so is / — ag for some a > small enough. Since 
the Ising model has the FKG property (see tt((/ — ctg)g) > ir(f — ag)wg so 
tK/s) — a7r (9 2 ) ^ ( 7r /)( 7r 5) — c^i^g) 2 so — airg 2 > 0, contradiction. 
Lemma|31guarantees there is an increasing g corresponding to 112', hence / must 
also correspond to fi2- 



3 Proof of the Main Result 

From now on we assume our graph G is a cycle, with vertices V — {1, 2, . . . , n} 
and edges [1, 2], [2, 3], . . . , [n,l]. All notation will be mod n, so n + 1 = 1. 
The coupling constants are J^+i = J; > 0. As usual, every vertex i £ V 
has a spin <7j G {+1,-1} and the probability of a configuration is given by 
7t((t) = Z^ 1 exp(^" =1 Jicrj<x; + i). We assume no external field and no boundary 
conditions. 

For technical reasons we assume Ji > 0; the case when some of the J's are 
zero can be easily obtained by limiting arguments. In particular, most results 
in this paper hold in the case when the graph G is a straight line (J n ,i = 0). 

We have already introduced increasing functions. Another important class 
of functions are linear functions, that is functions of the form f(a) — Yl^—i qiOi 
for some real constants qi. On the cycle we have 
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Lemma 5 If f is linear then Af is linear. 

Proof. We will need the standard hyperbolic functions sinh(a;) = l/2(exp(x) — 
cxp(— x)), cosh(x) = l/2(exp(x) + exp(— x)), tanh(x) = sinh(x)/ cosh(x). Since 
A is a linear operator, it is enough to consider f(a) = Oi. Then f(a J ) = Oi for 
j ^ i, and f(a l ) = -a u so 

(Af)(a) = A{a,a)f{a)+Y,A{a,ai)f{^) 



= (1 


3 

-^A{a,a j ))ai^ 




= 0"i 
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= (1 


- \jn)a l + (l/n)cr, 


t 1 /™) ^ . f i\ a * 

TT((J) + 7r(<7 1 ) 


= (1 


- l/n)a l + (l/n)<7j 


7r((7 l ) — 7r((T) 
7r((T l ) + 7r((T) 


= (1 


- l/n)(Tj + (l/n)o-j 


tanh(J i _i(TiO' i _i + Jj(TjC7, + i) 



Since tanh is an odd function, the second term does not depend on Oi , so 

(Af)(a) = (1 - l/n)<r t + (1/n) tanh(J i _i<r i _i + J;+i<7; + i) (4) 
and the tanh is easily expressed as a linear function 

tanh(J i _icr i _i + J i+1 a i+ i) = a<7j_i + ba i+1 (5) 

since Oi-\ and <Ji + \ only take the values +1 and —1, and both sides of the 
equality are odd functions. So we only need to check that a + b = tanh(Ji_i + 
Ji), a — b = tanh(Ji_i — Ji). After a bit of algebra, this yields 

a = sinh(2Ji_i)/(cosh(2 Ji_i) + cosh(2Jj)) 
b = sinh(2J 4 )/(cosh(2J,_i) + cosh(2J,)) 

and the lemma is proven. 

Hence we can look at the restriction of A to the space of linear functions. 
With respect to the basis {ci, a<i, . . . , <r n }, this is described by an n x n matrix 
L whose entries are all zero, except on the diagonal, and immediately above and 
below it. Let Sj = sinh(2Ji) and Cj = cosh(2Ji). Then we have 

L(i, i) = 1 — 1/n 
L(i,i-1) = {l/n)si-i/(ci-i+Ci) 
L(i,i+1) = (l/n)si/(ci-i + Ci) 
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Note that an eigenvalue of L is also an eigenvalue of A, so all must be real. Since 
all entries on the diagonal of L are equal, we can write L = (1 — l/n)I+ (l/n)M, 
where the diagonal of M is all zeroes, and the only non-zero entries are 

M(i,i-l) = Si- 1 /(*-i + Ci), M{i,i + l) = Si/(a-i + Ci) (6) 

Let the eigenvalues of JW be Ai > ... > A n . Then the eigenvalues of L are 
1 — (1 — \i)/n for 1 < i < n, and M and L have the same eigenvectors. It is 
easy to see M is irreducible. Hence by Perron-Frobenius, Ai > 0, Ai is a simple 
eigenvalue, and there exists a dominant eigenvector x (that is, corresponding to 
Ai) which has non-negative entries. In fact, the entries must be strictly positive, 
since we have Mx = Aix, so 

M(i, i - + M(i, i + l)x i+ i = XiXi (7) 

so if Xi = for some i, then also Xi+\ = 0, so x = 0, contradiction. But then 
f( a ) = ^2 x i (T i is a strictly increasing function of c, and also an cigcnfunction 
of A. Lemma 0] then implies that / corresponds to the second eigenvalue. We 
have obtained 

Corollary 6 ji% = 1 — (1 — Ai)/n, and /12 has a linear eigenf unction. 

Hence the relaxation time of the Glauber dynamics can be computed as 
T2 = l/(n(l— /12)) = 1/(1 — Ai), so from now on we can focus on the small matrix 
M. We will need an inequality which roughly says that its lead eigenvectors are 
not too far from a constant vector. 

Lemma 7 If x is a non-negative dominant eigenvector of M , then 

(ci - l)/si < Xi/x i+ i < (a + l)/si (8) 

Proof. The method could be called "cyclic induction." Note that cf — sf = 1, 
so (a — I)/ Si = Si/(a + 1) < 1. We have Mx = Aix, so 

Ai(ci_i -I- Ci)x t = Sj_i^_i + s^i+i Vi (9) 

Assume the left-hand inequality of JSJ holds for i — 1 so (c^_i — l)/si_i < 
Xi-x/xi. Combining this with ©, we get (since Ai < 1) 

(Ci_i + Ci)Xi > Ai(Ci_i + CijXi > (Cj_l - l)Xi + SiX i+1 

so (ci + l)xi > SiXi+i soxi/xi + i > Si/(ci + l) = (cj — l)/si. Hence the inequality 
holds for i as well; hence it holds either for all i or for none. But in the latter case 
we would get Xi/xi + i < (ci — l)/si < 1 for all i, hence x% < xi < . . . < x n < xi, 
which is impossible. 

The right-hand inequality can be proven in a similar way: if the inequality 
(cj + l)/si < Xi/xi+i holds for i — 1, it also holds for i; it cannot hold for all i, 
so it cannot hold for any i. Alternatively, the right-hand side can be deduced 
from the left-hand side by symmetry. 

An immediate consequence of this is 
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Proposition 8 //x is a non-negative dominant eigenvector of M , then Xi/xi+\ + 

Proof. The function fit) = t + 1/t is decreasing for t < 1 and increasing for 
t > 1, so it follows from Lemma [7| that f(xi/x i+1 ) < f((ci + l)/si) if £i > a^+i, 
and f(xi/xi+i) < f({ci — l)/si) if X{ < Xi+%. Both upper bounds are equal to 

To apply the last result we need an expression for the largest eigenvalue 
in terms of its eigenvector. Note that the product of the elements under 
the diagonal of M is equal to the product of the ones above the diagonal, 
so M is reversible in the following sense. Define q(i) = Ci + c,+i. Then 
q(i)M(i,j) — q(j)M(j, i)Vi,j. So if Q is the diagonal matrix with Q(i, i) = q(i) 
then Q 1 / 2 MQ~ 1 ^ 2 is symmetric, and has the same eigenvalues as M. Then by 
the usual extremal characterization, 

(Q^MQ-'/iv^v) 
Ai = sup 

(v,v) 

where the sup can be taken over all nonzero vectors v, and is attained at domi- 
nant eigenvectors of Q X I 2 MQ~ X I 2 . The numerator expands as 

(Q 1/2 MQ^/ 2 v,v) = 2ViV i+1 ^M{i, i + l)M(i + 1, i) 

= 2v l v i+1 s l / \J (cj_i + Cjjjcj + c i+ i) 

so if we finally substitute vi = (ci + Q+i) 1 / 2 ^ we obtain 

Lemma 9 Ai = sup y (Yl tyiUi+iSi) / (J2(yf + J/f+i)ci) and the sup is attained 
at dominant eigenvectors of M. 

Now let h(y, 3) = 2y l y i+1 s l ) / [Y^ivf + yf+i)ci) and let x = x(J) be a non- 
negative dominant eigenvector of M, chosen so x(-) is a smooth function in a 
neighborhood of J. This can be done because Ai is a simple eigenvalue; see, for 
example, 0, Chapter 2.1. We have Ai(J) = /i(x(J), J). Because for fixed J, 
sup/i(y, J) is attained at y = x(J), we have 

dh/dy k (x(3),3) =0 Vfc (10) 

so by the chain rule 

OXt/OJi = dh/dJ, + J2(dh/dy k ){dx k /dJi) = dh/dJi (11) 
fc 

Hence 

d\ 1 /dJ i = dh/dJi(x,3) 

= 2(2xiX i+1 Ci - (x 2 + x 2 +1 )sih(x,3))/^(x 2 + x 2 +1 )ci (12) 
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From Proposition|Sl (x^+xl +1 )/{2xiX i+ i) < Ci/si and since h(x, J) = Ai < 1 
it follows that dX\/dJi > 0. This completes the proof of Theorem 

4 Further Inequalities 

We can get more precise bounds on the derivatives of Ai. Indeed, equation 1|I2[I 
and Proposition IS] yield easily 

dX^d.h > 2(1 - X 1 )(2x t x z+1 c z )/J2(^+^+i)c t (13) 

hence 

^2(s i /c i )dX 1 /dJ i > 2(1 - Ai) ^2(2xiX i+1 Si)/ y^(x- + xf +1 )a 

= 2(1-A 1 )A 1 (14) 

So we can obtain an interesting inequality by differentiating along a par- 
ticular direction. To do that, consider functions Ji(t) with J^(0) = Ji and 
J'i{t) = tanh(2 Ji(t)). This differential equation is easy to solve and yields 

Ji{t) = (1/2) sinh _1 (exp(2t) sinh(2Jj)) (15) 

(here sinh -1 denotes the inverse hyperbolic function). Then the function X(t) = 
Ai(Ji(<), . . • , J n (t)) satisfies A' > 2A(1 — A). This becomes even simpler if we 
substitute 

t* = t 2 - 1 = A/(l - A) (16) 

We obtain r*' > 2r* so (logr*)' > 2 and r*(t) > exp(2i). Hence, with a 
slight abuse of notation (by /(J) we mean that / is applied to all entries of J), 
we have 

Proposition 10 r*((l/2) sinh~ 1 (exp(2t) sinh(2J))) > exp(2t)r*(J)VJ, t > 0. 

For large x, sinh _1 (a;) is close to log(a;), so we expect Ji(t) to be approxi- 
mately Ji + t. In fact it is easy to prove that exp(2£) sinh(2u) < sinh(2t + 2u), 
so Ji(t) < Ji +t, so from Theorem \I\we get 

Theorem 11 t*(J + 1) > exp(2t)r*(J) for all t > 0. 

As above, J-M means t is added to all entries of J. Numerical computations 
confirm that the factor 2 in the exponent is optimal. It would be nice to have 
a similar inequality where the elements of J are multiplied by t, so that we can 
estimate the effect of a change of temperature on the relaxation time. Note, 
however, that t J > t + 3 for large t, so the inequality we have proven is in some 
sense stronger. 
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5 Low Temperature Asymptotics 



Introduce the inverse temperature /3; the coupling constants are now j3J\ , . . . , (3J r 
We are interested in T2(/3J) for large f3 and fixed J. This turns out to be asymp- 
totically exponential in (3. More precisely, let J max — max Ji and J m in — min Ji. 
We have 

Proposition 12 The following holds as (3 — > oo . ■ 

limr 2 (/?J)/cxp(C / 9) = D 

where C = 2(J max + J mm ) and D = (l/2)#{z : J t = J ma x}/#{i ■ Ji = Jmin}- 

Proof. As before, let x be a dominant eigenvector. Lemma [7| guarantees that 
the ratios Xi/xi+± are close to 1 for large j3. That is, given e > 0, 1 — e < 
xi/xi + i < 1 + e for (3 large enough. We can set x\ = 1 so we get (1 — e) n < 
Xi < (l + e) n Vi. 

If we apply Lemma|21for y = 1 we obtain Ai > y]sj/y]cj so 

r 2 = l/(l-A 1 )>^c i /^(c i - Si ) (17) 

If we apply Lemma [5] for y = x and substitute t-i = 1/(1 — Ai), we obtain 

t 2 = y^Cj(xj + a;f +1 )/[^Ci(gj - x i+ i) 2 + ' s ^2x i x i+1 (c i - Sj)] 

< ((l + e )/(l- e )) 2 «^ Cl /^( Cl - Sl ) (18) 

The ratio X) c i/XX c i ~ s *) ^ s easily shown to be asymptotically Dexp(C/3). 
Since e is arbitrary, the proof is complete. 

Acknowledgments. Many thanks to Yuval Peres, my thesis advisor. His 
guidance was instrumental in obtaining these results. I also thank him for per- 
mission to include Lemmas |3 and 0| Thanks also to Noam Berger, Elchanan 
Mossel and Dror Weitz for their comments. 



References 

[1] Fortuin, C. M., Kasteleyn, P. W. and Ginibre, J. (1971) Correlation in- 
equalities in some partially ordered sets, Comm. Math. Phys 22, 89-103. 

[2] Kato, T. (1984) Perturbation Theory for Linear Operators, Springer- 
Verlag. 

[3] Kenyon, C, Mossel, E. and Peres, Y. (2001) Glauber Dynamics on Trees 
and Hyperbolic Graphs, Forty-Second IEEE Symposium on Foundations 
of Computer Science (FOCS). 



9 



[4] Liggett, T. M. (1985) Interacting particle systems, Springer- Verlag. 

[5] Martinelli, F. (1998) Lectures on Glauber dynamics for discrete spin mod- 
els. Lectures on probability theory and statistics (Saint-Flour, 1997) 93-191, 
Lecture Notes in Math. 1717, Springer, Berlin. 



10 



